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We present examples of non-Gaussian statistics that can induce bispectra matching local and 
non-local (including equilateral) templates in biased sub-volumes. We find cases where the biasing 
from coupling to long wavelength modes affects only the power spectrum, only the bispectrum or 
both. Our results suggest that ruling out multi-field scenarios is quite difficult: some measurements 
of the scalar correlations, including the shape of the bispectrum, can be consistent with single¬ 
clock predictions even when cosmic variance from super-horizon modes is at work. Furthermore, if 
observations of the density perturbations rule out single-clock inflation, we will face a serious cosmic 
variance obstacle in drawing any further conclusions about the particle physics origin of the scalar 
fluctuations. 

PACS numbers: 


I. INTRODUCTION 

The primordial curvature perturbations are our primary source of information about the inflationary era of the 
universe. Ideally, we would like to extract an understanding of the particle physics responsible for inflation from the 
details of the correlations observed in the Cosmic Microwave Background (CMB) inhomogeneities and the large scale 
structure. However, we observe only a finite volume of the universe and have no reason to expect that the current 
size of that volume is anything special: there are likely to be at least some super horizon modes with more or less 
the same properties as those that have already re-entered the horizon. That supposition is critically important in 
comparing statistical observations to theory. In the case of exactly Gaussian fluctuations the finite size of our Hubble 
volume is the origin of the familiar cosmic variance uncertainties in the power spectrum fit. 

In non-Gaussian scenarios that couple Fourier modes of very different wavelengths there is additional cosmic variance 
from the possibility that all our observations are biased compared to the mean predictions of some inflationary theory. 
The bias from super horizon modes is completely unmeasurable, but the qualitative conclusions we draw about the 
origin of fluctuations can change when we allow for it pm Fortunately, no such mode coupling has been detected 
within our universe yet, but current observational limits (e.g., from the Planck satellite results m and the most 
recent constraints from quasars m) do not rule out this possibility. Indeed, many inflation scenarios predict some 
degree of non-Gaussianity and the details of the correlations would ideally provide a means to distinguish qualitatively 
different primordial physics. In weighing observational evidence for or against any inflation model that couples modes 
of different wavelengths we must include these “Super Cosmic Variance” (SCV) uncertainties. In addition, if any non- 
Gaussianity is detected, we must be sure we understand how to draw robust conclusions about the set of primordial 
universe models consistent with those measurements. This is a non-trivial task because it is not clear if there are 
types of non-Gaussianity inflation cannot generate and because it is far from obvious how the statistics in biased 
sub-volumes are related to those of parent distributions with arbitrary non-Gaussian fluctuations. 

Here we consider several scenarios motivated by predictions from currently studied inflation models and demonstrate 
how the presence of correlations beyond the bispectrum can alter the shape of lower order correlation functions in 
biased sub-volumes. We will demonstrate that although observations could prove that the source of the primordial 
fluctuations was not single-clock inflation, constraints on or detections of the shape of the bispectrum can be consistent 
with single-clock predictions even when super cosmic variance is at work. 
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A. The model 


An attractive discriminating feature of inflation scenarios is the behavior of the squeezed limit of the primordial 
bispectrum: that is, how significantly it couples two short wavelength modes to one long wavelength mode (represented 
by a squeezed triangle in momentum space, with side lengths k\ = ki <C fc 2 ~ £3 = h) m ■ This limit also indicates 
how significantly the bispectrum can cause the power spectrum in biased sub-volumes to differ from the global power 
spectrum (for a concrete example entirely within our universe, see the discussion of non-Gaussian halo bias • KlI 1141 1. 
A bispectrum of the type generated by single clock inflation, which primarily couples modes of the same wavelength, 
will give a negligible shift to the power spectrum regardless of the realization of the long wavelength modes. A local 
type coupling, on the other hand, can give an interesting amplitude shift HE]- 

Beyond the bispectrum, all higher order correlation functions can contribute to the biasing of lower order statistics 
when some modes have longer wavelengths than the size of the spatial region of interest. Scenarios where the non- 
Gaussian field is any local but non-linear function of a Gaussian field have the nice property that the bispectrum is 
of the standard local type (up to at most logarithmic corrections) even in biased sub-volumes. And, in sufficiently 
biased sub-volumes the observed statistics of any local, non-linear function (polynomial and without derivatives) of a 
Gaussian field will be those of a field with the local ansatz [2]. 

Here, we explore the question of biased statistics for scenarios with non-local bispectra. Our goal is to understand 
how observations of a necessarily limited set of correlation functions can constrain the space of models of the primordial 
universe. In particular, the single-clock inflation consistency relations H2 G3H22] indicate how one could rule out 
single clock inflation, but would a detection of £ al = 0 and 7 ^ 0 confirm single clock inflation? To address this 

question we consider bispectral shapes characterized by their degree of divergence with the long wavelength mode, 
hi, in the squeezed limit, including k ; _1 (equilateral type), fc ; -2 (sometimes called ‘orthogonal’ type 1 ), and k^ ’’ (local 
type). 

To consider non-Gaussian scenarios with the desired 71 -point correlations, we build the field $ from a series of 
nonlocal functionals of a Gaussian random field 4>(x), 

$[<K X )] = <K X ) + /NL$2 [</>( x )] + ffNL$3[0( x )] + • • • , (1) 

where the subscript on indicates how many copies of the Gaussian field appear in the term. The <f> ra will generate 
the connected parts of the tree-level correlations at order n + 1 and higher, and we require ($„) = 0. We assume the 
Gaussian field is homogeneous and isotropic. Its statistics are completely determined by the power spectrum, which 
we take to be scale-invariant for simplicity: 

<^(k 1 )^(k 2 )) = (27r) 3 <5 3 (k 1 +k 2 )P^(fc 1 ) (2) 

Pm = 2tt 2 ^ . (3) 

We will primarily work with the Fourier transform of <f>: 


<f>(k) — </>(k) + /NL < l ) 2(k) + «/NL^3(k) + ... 
/nl f d 3 p 1 


<M k ) + ^ J J d3 Pi [ < Mpi) ( KP2) - W>(Pi)</>(P2))] IV 2 (Pi,P2,k)<T 


(k-pi -p 2 ) 


9 NL 


3!(2tt) 6 


n /* 


pi 


</>(Pi)<MP 2)</>(P3) - < KPi)(0(Pj) < / , (Pfc)> 


i =1 
k^jy£i 


N 3 (Pi,P2,P3,k)<5 3 (k-^p f ) (4) 


1=1 


The kernels N n (jpi, p 2 ,..., p n , k) are symmetric in the first n entries. The subtracted terms inside the square brackets 
maintain (<f> n ) = 0 and ensure that only connected parts of the n-point functions are generated by each term. Loop 
corrections to the power spectrum from the non-Gaussian terms go like (/nlA^) 2 and (^nlA^,) 2 and the one-loop 
correction to the bispectrum goes like (?nlA 2 (relative to the C>(/nl) tree-level bispectrum). Weak non-Gaussianity 


1 A different, also very useful definition of ‘orthogonal’ comes from looking at possible scale-invariant shapes produced by single clock 
inflation m- With that definition, the orthogonal template diverges as 1 . 
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is defined by requiring those quantities to be <C 1, and if ^nl < ^(/nl) the 1°°P correction to the power spectrum 
from the three point function is parametrically larger than that from the four-point. 

To generate our full set of desired correlations, the 4>„ must be non-local. It is notationally easier to first generate 
the appropriate set of terms in real space, so we consider <E> n of the form 

= 9^"- 1 (5 /32 "- 2 0(...(5 /32 ^5 /3l 0)))(x) , (5) 


where the fti can be negative. To generate the corresponding Fourier transformed terms, we define derivative operators 
acting on </>(x) based on the corresponding momentum space behavior, 2 


a n </>(x) = 


d 3 k 


k n cj)( k) 


Ah-: 


( 7 ) 


A $ 2 that can generate the equilateral, orthogonal or local bispectral templates was derived by Scoccimarro et al in 
123) . In the next section we review the choice of quadratic terms, taking a slightly different perspective but arriving 
at the same $ 2 as [24]. In Section III we generalize the procedure for generating the quadratic terms and consider 


a d> 3 designed to induce only the terms in <J> 2 in biased sub-volumes. In Section [TV] we examine subsets of the cubic 
terms that restrict the bispectra induced in biased sub-volumes to local, equilateral, or orthogonal type and consider 
the implications for theory and observations. In Section [U] we conclude and comment on how our analysis might be 
further developed. 


II. QUADRATIC TERMS 

In this section we derive an expression for $ 2 using a procedure that can be easily generalized to the cubic term, 
<I> 3 , and beyond. Although our motivation and procedure are slightly different, the result reproduces the expressions 
derived by Scoccimarro et al [24] . 

We limit our considerations to bispectra that diverge in the squeezed limit as one, two, or three inverse powers of 
the long wavelength momentum. Since a local quadratic function, <fi 2 (x), yields the most divergent bispectrum, we 
need to include inverse derivative operators in order to reduce the divergence. We therefore start with the following 
family of terms quadratic in the Gaussian field 


d a3 (d a2 (t>d ai (l)), 


( 8 ) 


with restrictions on the af. 

• Yl a i = 0- This condition maintains scale-invariance. 

• oq i2 > 0. Together with the previous condition, this automatically sets <23 < 0. This condition ensures that 
the infrared (IR) sensitivity to <f>i is no stronger than that of local non-Gaussianity. In other words, this rule 
ensures that the squeezed limit of the bispectrum does not grow like k7 i or more for small /q. 

• | oti | < 2. This restriction generates the minimal set of terms required to produce equilateral (fc ; _1 ), orthogonal 
(fc; -2 ) or local (fc ; _3 ) type behavior (and corresponds to to setting u = s = 0 in [24]). However, there are 
certainly additional terms with \cti\ > 2 that will also generate bispectra with the same squeezed limit behavior. 
We comment on relaxing this condition in Section |V| 

The generic quadratic functional with terms obeying these constraints is 

4> 2 [<^>(cc)] = [ai<^> 2 + a 2 d~ 1 ((/)d(f>) + a 3 d~ 2 (4>d 2 </>) + a4<9 _2 (<9<(>) 2 ] - [E.V.] . (9) 


2 Note that d n is not in general a genuine derivative operator because it does not obey the Leibniz rule, 

d 3 p 


[d n (?ii(/>2)](k) = fc- / 

{<fiid n fo + 0 2 0"0i] (k) = |^0l(p)02(k-p)(p" + |k-p|’ 1 ) 
-A d n (<f)l(j)2) / (pld n (f)2 + 

Our notation is related to the one used in Scoccimarro et al m by: (<9 2 )here = (~V 2 ) there• 


(6) 
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where — [E.V.] indicates that the expectation values of the terms should be subtracted. The corresponding JV 2 kernel 
(which in this case depends only on the magnitudes of the momenta) is 


N 2 (pi,P 2 ,k) = 2a 1 +a 2 


P 1 +P 2 , p\+pI . „ P 1 P 2 


A generic homogeneous and isotropic bispectrum for the potential $ can be written as 

($(k 1 )$(k 2 )$(k 3 )) = (27r) 3 A 3 (k! + k 2 + k 3 ) B(k u k 2 , k 3 ) ■ 
Our ansatz gives a bispectrum 

B${k 1 , fc 2 , fc 3 ) = /NL-P 0 (fci)f^(fc 2 )A r 2 (fci, k 2 , k 3 ) + eye. 


( 10 ) 


( 11 ) 


( 12 ) 


where there are 2 additional cyclic permutations. Notice that both the a 2 and a 3 terms in Eq.(10) generate terms of 
the type P(fci) 1 ,/ 3 P(fc 2 ) 2 / 3 P(fc 3 ) in the bispectrum. 

We have already factored out an overall amplitude, /nl, from $ 2 . With the usual convention / NL = 
6 P$(fc, k, fc)/Pj(fc), this leaves a normalization condition for the coefficients of the individual quadratic terms: 


Cbl + CL 2 + U 3 + 04 — 1 • (13) 

The constraints on the a* so far allow us to restrict consideration to a subset of possible non-Gaussian fields based 
on the behavior of the bispectrum. However, any quadratic term will also contribute to the power spectrum of the 
full, non-Gaussian field 


($(k 1 )$(k 2 )) = MkiMka)) + / 2 L (4> 2 (k 1 )4> 2 (k 2 )) + ... (14) 

where the dots contain contributions from <f > 3 and higher. The contributions to P* (fc) from the terms in <f > 2 contain 
an extra integral over momenta p and go like [N 2 (p, |k — p| , k)] 2 . The contributions from the a 3 and 04 terms go like 
1/fc 4 times a divergent integral over momenta p (this is easiest to see by considering iCp). We can force this badly 
behaved contribution to the loop to vanish by setting 


04 = —a 3 . 


(15) 


We will impose this condition from now on to remove 04 from all expressions, 
the choice of the coefficient t for the equilateral and orthogonal cases in [23j . 


Insisting on Eq.(15) is equivalent to 


A. Recovering the standard bispectral templates 


Frequently used bispectral templates with fixed degree of divergence with the long wavelength mode are 

Plocal = 2/)? L Cal (P 0 (fc 1 )P^(fc 2 )+P ?i (fc 1 )P 0 (fc 3 )+P^(fc 2 )P 0 (fc 3 )) (16) 

Pequil = 6 /^ q L uil [-P 0 (A: 1 )P 0 (fc 2 ) + 2 perm. - 2 (P 0 (A: 1 )P 0 (fc 2 )P 0 (fc 3 )) 2/3 
+P 0 (fci) 1 / 3 P 0 (fc 2 ) 2 / 3 P 0 (fc 3 ) + 5 perm.] 

Perth = 6/£ r L th [-3P 0 (fci)P 0 (fc 2 ) + 2 perm. - 8 (P 0 (fci)P 0 (fc 2 )P^(fc 3 )) 2/3 
+ 3 P 0 (fci) 1 / 3 P^(fc 2 ) 2 / 3 P 0 (fc 3 ) + 5 perm.] 


The best constraints on the amplitudes /nl of these templates come from the Planck satellite: m, which limit 
/i°cai = 2.7 ± 5 . 8 , / aqail = —42 ± 75, and /§ r L th = -25 ± 39 at the 68.3% confidence level. 

The ansatz in Eq.® clearly contains the local ansatz 


Local Bispectrum : a 2 = a 3 = 0,4 = 0 . 


(17) 


To see the conditions that recover the orthogonal and equilateral templates, insert the N 2 kernel from Eq. (10 1 into 
the general expression for the bispectrum, Eq. (12), and take the squeezed limit k\ = ki fc 2 ,fc 3 = k s : 


lim B(k s ,k s , h) = fNi^P^k^P^ks 


k t ^k 


(4ai + 2a 2 + 2a 3 ) + (2a 2 — 4a 3 ) + (2ft2 + 2a 3 ) ''j + 2ai 


k ' 


■ki¬ 


rn 
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where we have already used <24 = —< 23 . The equilateral bispectrum scaling k l 1 requires that the coefficients of both 
the fc ; -3 and the k~[ 2 contributions vanish. Since P{k{) oc kf 3 , these conditions are, respectively 


4a 1 + 2a 2 + 2 a 3 = 0 and 2 a 2 — 4«3 = 0 . 


(19) 


Including the normalization condition (13), these are sufficient equations to uniquely fix the ai 
recovers the equilateral template: 


The result exactly 


Equilateral Bispectrum : ai = —3, a 2 = 4, 03 = 2, a 4 = — 2 . (20) 

The orthogonal bispectrum scaling k7 2 requires only the coefficient of k ; -3 to vanish, which is not sufficient in¬ 
formation to completely fix all the ai. That is reasonable since the orthogonal shape was originally defined not by 
its squeezed limit but by minimizing its overlap with the other templates over a range of momentum configurations. 
To generate the standard orthogonal template we can add an additional condition, comparable to the orthogonality 
condition required by [23], 

a° rth = 3 a® quil (21) 

so that all the are fixed: 


Orthogonal Bispectrum : ai = —9, a 2 = 10, 03 = 8 , 04 = —8 . 


( 22 ) 


In the next section we will introduce the long-short wavelength split and see how the conditions in Eq.(19) that fix 
the squeezed limit behavior of the bispectrum can also be directly read off of the real-space expression in Eq. 


B. The long-short wavelength split 


We are interested in the effect of long wavelength background modes on the statistics measured in spatial sub¬ 
volumes. To that end, we split the field $ in Fourier space at a scale &*, 

$ s (k) = $(k)e(fc-fc 0 

$*(k) = $(k )©(fc*-fc), (23) 


where 0 is the step function. This is only an approximation to splitting the field with a top hat in real space, but 
our results will not depend on this distinction. Also, in what follows, we will mainly consider momenta far away from 
the scale fc*, and ignore complications arising from those close to fc*. Applying the 0-function on the right hand side 
of $(k), defined in Eq. © , gives 

$ s (k) = ^(k) + / NL $^(k) + ... (24) 

$*(k) = ^(k) + /NL^ ) (k) + ... (25) 


where 0 s (k) and <pi( k) are defined analogously to Eq. (23). Writing out the effect of the 0-function on the quadratic 
term gives, for the short wavelength quadratic piece, 


4 s \k) 


9 j ( 9.^)3 J d3pi d 3 P2 N 2 {Pi,P2^)[^(Pi)^(P2) - [E.V.]] 5 3 (k- Pi - p 2 ) 0 (|pi +p 2 | - fc«) 

(27i')3 / d 3 pid 3 p2 A^ 2 (pi , P 2 , k)[<^> s (pi)(^> s (p 2 ) T 2<j6 s (pi )</ , z (p 2 ) [E-V.]] <5 3 (k pi p 2 ) (26) 

|pi+P2|>k» 


Here, the second line is obtained from by considering the separate momentum regimes for pi and p 2 that satisfy the 
step function: 0(|pi+p 2 |-fc*) = [0(pi-fc*)0(p 2 -/c») + 0(pi -fc*)0(fc*-p 2 ) + 0(fc*-pi)0(p 2 -fc*)]0(|pi + p 2 | -fc»). 
In the same manner we can write $ 2 ^ (with some redundancy) as 

$2 } ( k ) = 9 |( 9^.) 3 J d 3 pi d 3 p 2 N 2 (pi,p 2 , k)[0 s (pi)</> s (p 2 ) + cj)i(pi)(j)i(Tp 2 ) — [E.V.]] 5 3 (k — pi — p 2 ) (27) 

|pi+P2l<k, 


Notice that in the first term of Eq.(27), the angle between short scale modes pi and p 2 must be nearly 7 r so they sum 
to a long mode with magnitude k below the splitting scale. 
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C. The field observed in biased sub-volumes 

We now examine the difference in the global and local statistics of the power spectrum when the quadratic terms 
above are included. Very generally, the observed field, restricted to points x within a sub-volume Vs, will take the 
form 


$° bs (x)U 6 y s = (</> + $? CV ) + /nl($2 + V ) + <7Nl($3 + V ) ■ • • , (28) 

where the contributions from the Super Cosmic Variance (SCV) terms, < f > ® cv , vanish in sub-volumes where all long 
wavelength inodes take their mean value (0). Since an observer with access to only a single sub-volume cannot separate 
the mean contributions from the SCV contributions it is more natural to write the right hand side of the equation 
above in terms of the observer’s linear field, y(x), and kernels, and amplitudes: 

$obs L\( x )]UeVs = x( x ) + /nl'MxM] + <7 NL i> 3 [x(x)].. • , (29) 


Note that the standard way of defining the normalizations /nl and on ly works for scale-invariant contributions 
from the SCV terms. 

For the set of quadratic terms we are considering here, using the kernel N 2 from Eq.(10) in the expression for the 
short wavelength field, Eq. (26) gives 

$ s (k) ~ </> s (k)j 1 + / NL (2ai + a 2 + a 3 )/>; + (a 2 - 2 a 3 )dcf>i^ + a 3 d 2 (f>i 


+ 


/nl 

2!(2tt) 3 


d 3 p 1 d 3 p 2 N 2 (p 1 ,p 2 ,k) 


Pi,P2>k, 


4*s (pi)^s (P 2 ) - [E.v.] 


<5 3 (k pi p 2 ) 


(30) 


where the approximation comes because we have taken k —p 2 ~ k in the </> s (pi)</>;(p 2 ) term in Eq.(26). The functions 
of the long wavelength modes are 


<t>l 


> 


<Kp), d(j)i = 


d 3 p 
(2t rp 


P^{ P), d (j)i = 


d 3 p 

(2tt) 3 


P 2 H P) 


(31) 


where A is an infrared scale corresponding to the longest wave mode that exited the horizon during inflation. We 
will treat the quantities in Eq. (31) as constants in any particular sub-volume. While this is not exactly true, since a 
nonzero gradient difii implies non-constant />;(x), we show below that this difference is small compared to the effect 
of the average value of </>j. 

The observed linear term in a biased sub-volume is shifted from the original </>(fc) by a term whose amplitude and 
scale-dependence depends on the bispectrum and the bias: 


Xc(k) = 0 s (k)[l +/ NL $f cv (fc)] 

= </> s (k)jl + / NL (2ai +a 2 +a 3 )(/)i + (a 2 - 2a 3 )<9</>; ^ +a 3 <9 2 </>;^j 

Then, using the observed linear field everywhere, the observed non-Gaussian field is 
$ ohs (k) = X g( k) 


(32) 


(33) 


/nl 

"2!(2tt) 


d 3 pi d 3 p 2 


V 2 (pi,p 2 ,k) 


[1 +/nl^ CV (pi)][1 +/nl^ CV (p 2 )] 


[xg(Pi)xg(P 2) - [E.V.]]d> 3 (k - pi - p 2 ) . 


The power spectrum and bispectrum observed in sub-volumes can be computed as usual from this expression, and 
will differ in amplitude and scale-dependence (shape) from the corresponding quantities in unbiased sub-volumes. 
The observed power spectrum is shifted from the input power spectrum by (assuming weak non-Gaussianity) 


P' 


iobs 


P° bb (k) — P(j,{k)\ 1 + /nl (2ai + a 2 + a 3 )</>z + (a 2 — 2 a 3 )d(f>i— + a 3 d 2 (j)i—,^ 


(34) 


This effect was discussed and quantified in detail for the case of local non-Gaussianity in [4]. There it was shown that 
the correction proportional to (f>i can be large and interesting. The other terms here, proportional to d(f>i and d 2 (f>i are 
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sensitive to only a very small range of modes beyond the horizon. If the universe was inflating for N e e-folds before 
the mode of scale k * exited the horizon, the super cosmic variance contributions from non-local bispectra are of order 


mi) 2 ) 

k* 

((d 2 m) 

K 


1 

fcf 


1 

kf 




(35) 


Here we have assumed for simplicity that the super horizon power spectrum is = Ao(fc//c*) ns 1 with constant 
spectral index n s and used N e = ln(k*/ A). Notice that with fc* = Hq the second quantity in Eq. ([35]) is of the order 
of the contribution of fluctuations to the sp atial curvature, (f 2 |). 

Figure [T] shows that the quantities in Eq.(35) are very small, of order 10 -5 . As expected, there is no appreciable 
cosmic variance shift to the power spectrum, Eq. (34), from non-local bispectra. Figure IT] assumes that the power 
spectrum is consistent with the current best fit from the Planck satellite (Planck- 1 -WP) [25] (A %(k) = 3.98772 x 
10^(fe/fc*)"*" 1 , n s = 0.9619) and uses fc* = 0.05Mpc -1 as the reference infrared scale. ( H 0 is the appropriate scale 
for considering our current universe.) The strong suppression with extra powers of k in the integrands makes these 
quantities insensitive to small variations in the spectral index. Specifically, there is no considerable difference between 
flat spectral index (n s = 1) and n s = 0.9619 (in contrast to the significant enhancement (c ^) 1 / 2 receives from a red 
tilt over sufficiently many super horizon e-folds). 



FIG. 1: Graph of the RMS amplitude of long-wavelength super cosmic variance terms from non-local bispectra, -j)~ {(drfn) 2 } 1 ^ 2 
and pj((d 2 ^) 2 ) 1 / 2 . This plot assumes that the power spectrum on scales larger than k » = 0.05Mpc -1 remains consistent with 

Planck satellite observations: A| (k) = 3.98772 x 10 -9 (fc/fc*) n5_1 with n s = 0.9619 [25]. Notice that the plotted quantities 
are not sensitive to more than about an e-fold of inflation and remain C2(10~ 5 ) so, as expected, only the local bispectrum can 
generate a significant super cosmic variance contribution to the power spectrum. 

Although these terms have a substantial scale-dependence, the maximum shift to the spectral index at the CMB 
pivot point is of order 1(E 4 , which is within the current 68 % confidence interval from measurements by the Planck 
satellite. (Notice that allowing mild scale-dependence in a 1; 02 or a 3 would generate interesting shifts to the spectral 
index. The example of scale-dependent ai was discussed in [5].) 
















The bispectrum observed in biased sub-volumes is also shifted. For the exact local ansatz, the change in the observed 
level of non-Gaussianity was discussed in mm- When the super cosmic variance shifts to the power spectrum are 
scale-dependent, the shape of the observed bispectrum changes: 


B ohs {k 1 ,k 2 ,k 3 ) 


P x (fci)P x (fc 2 ) 


/nl^W(^i, k 2: k 3 ) 

[1 + /NL$? CV (fci)][l + /nl<1>? CV (A)] 


+ eye. 


(36) 


However, as with the scale-dependent corrections to the power spectrum, the changes in shape are too small to be 
observationally relevant. We will find much more interesting effects from the inclusion of cubic terms in the ansatz 
for the non-Gaussian field, which we turn to next. 


III. CUBIC TERMS 


In this section we establish a method to generate cubic terms that induce the four quadratic terms of Eq. © in 
biased sub-volumes. 


A. Generation of cubic terms 

In general, we consider cubic terms that take the form 

d^ {d pi (j)d p3 (, d 02 ^ </>)). (37) 

Similarly to the rules we used to limit the quadratic terms, we consider only cubic terms that satisfy the following 
restrictions: 

1 . Y A = 0 to maintain scale-invariance, 

2 . P 1,2,4 > 0 to ensure that the induced quadratic terms do not depend on d~ 1 (f>i or higher inverse derivatives, 
controlling sensitivity to the infrared (IR) scale. 

3- | A | < 2 restricts to a minimal set of terms to consider. Note that conditions 1-3 imply (3 5 < 0. 

4. /3-| + B 2 + p 3 > 0 to ensure that the induced linear terms are no more sensitive to the IR scale than in local 
non-Gaussianity (logarithmic sensitivity). 

5. At least one of A, / 3 2 , or A must be zero and the {/A, A} pairs should be drawn from the same set of values 
as the {ai,a 2 } pairs in Eq.([9|. This ensures that all the leading order quadratic terms induced in biased 
sub-volumes are one of those in Eq.([9]). 

As a guide to understanding the origin of these restrictions, an example of the quadratic terms induced in biased 
sub-volumes by a cubic term is given in Appendix [B] There are 18 cubic terms, listed in Appendix [Aj which satisfy 
these five rules. As with the quadratic terms, it is certainly possible to consider a larger set of terms that have 
interesting behavior in biased sub-volumes. In Section [IV| we discuss the consequences of the minimal set of terms and 
also comment on the role of some terms we have discarded. Appendix |A| also gives the cubic kernel, iV 3 (pi, p 2 , P 3 , k), 
from these 18 terms. 

The addition of this cubic functional adds a leading order trispectrum to the model: 

(4>(ki)$(k 2 )$(k 3 )$(k4)) = (2?r) 3 (5|)(ki + k 2 + k 3 + k 4 ) T$(ki, k 2 , k 3 , k 4 ) 

T$(ki,k 2 ,k 3 ,k 4 ) = 5 N L-P 0 (fci)R^(A)R 0 (fc 3 )A r 3 (ki,k 2 ,k 3 ,k 4 ) + eye. (38) 

where there are three additional terms in the -l-cyc. Notice that our definition of ^nl as the amplitude of the 
trispectrum goes beyond the usual definition. The typical conventions for defining the amplitudes of trispectra in 
various momentum configurations are mm 

Standard trispectra conventions : 5 A andard = 1 lim k3 ’^ x 

WNL 6 ki-K) 

= 1_9_ r$(ki,k 2 ,k 3 ,k 4 ) 

TNL “ 4 25 Iki+kaj-^o P$(|k 4 + k 2 |)P$(fc 1 )P$(fc 3 ) 1 ’ 

Because we have imposed |A| < 2 (and specifically /3 3 > —2), all of our trispectra have tnl = 0. Current constraints 
on the trispectrum from Planck satellite data are (allowing both standard shapes to be non-zero) tnl = 0.3±0.9 x 10 4 , 
ffNL ndard = ~ 1-2 ± 2.8 x 10 5 H] 68 % CL, or assuming only one non-zero template at a time r NL < 2800 at 95% 
CL m, £$f mdard = -1.3 ± 1.8 x 10 5 at 68 % CL [2S]. 
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B. Constraints on cubic terms from their contributions to the power spectrum 


As with the quadratic terms, we impose additional restrictions based on the behavior of (classical) loop corrections 
from the cubic terms, requiring that UV divergences are not stronger than the log divergence of the tree-level case. 
The cubic field $3 contributes to the power spectrum at order (from (<& 3 (ki)$ 3 (k 2 ))): 


<SP$ 3 (fc) — l 


r d 3 pi d 3 p 2 
J (2tt)3 (2tt)3 


d 3 P3 N$(p 1 ,p2,P3,k)P(p 1 )P{p2)P(p 3 ) 5 3 (k - Pi - p 2 


P 3 ) ■ 


(40) 


Since this contribution has two momentum integrals (two classical loops), there are several different momentum 
configurations to consider. The resulting constraints (details appear in Appendix Jc| are 


bio = b 14 = bio = bis = 0 

b 3 + b± = 0 
bo b? — 0 
bi2 + bi3 = 0 
bn + bn + bn = 0 

2(65 + bg) + bo + bg + bi2 = 0 


(41) 


After imposing these restrictions, we are left with a 9 parameter set of cubic terms to explore. Corrections to the 
bispectrum from the cubic term (proportional to /nl<?nl at lowest order) do not give any additional constraints that 
are not covered by conditions imposed on the quadratic term. 


C. The field observed in biased sub-volumes 


To work out the statistics observed in biased sub-volumes of a Held described by Eq. 0 and the cubic functional 
from the previous section, we again split 4>(k) into long and short wavelength components: 


$a(k) = <Mk) +f NL ^ s) (k) + g NL ^ s) (k) (42) 

$z( k) = Mk) + fNL^k) + g NL ^ 1 ) (k) (43) 


where 


sO/'M - A 
3! 


^3 (k) = 


d 3 Pi d 3 p 2 3 


(27t) 3 (27t) 3 


d P3 ^ 3 (Pi, P 2 , P3, k) <5 J (k - Pi - p 2 p 3 ) B(|pi + P 2 + p 3 | - K) 


0(Pl)</>(P2)0(P3) - ^(P0(</ > (Pj)</ , (Pfc)) 


*=1 


+ eye 


(44) 


^ ){k) = h 


d 3 pi d 3 p 2 3 


(27t) 3 (27t) 3 


d P3^3(Pl,P2,P3,k) <T(k - pi - p 2 - p 3 )0(fc>. - |pi +P 2 +P 3 |) 


0(Pl)<^(P2)</>(P 3 ) - <KPi)(<KPj)0(Pfc)) 


i= 1 


eye 


(45) 


The step function in Eq.(44) allows several different types of terms (analogous to the two terms in Eq.(26)). In 


particular, there is a term with all three fields (j> s which ensures that all sub-volumes have the same cubic term as the 
parent volume. There is also a term with a single <pi that contributes to the quadratic term in biased sub-volumes. 
In other words, the quadratic term that describes the non-Gaussian field in the sub-volume (see Eq. (28)) is 


(/NL$2)° bS = /NL(d>2 + ®I CV ) 


(46) 
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where 


/nl^I Cy = 5 nl 


(tn{ ( 3 &i + b 2 + 63 + b 8 + bn) 4 > 2 s + (& 2 + b 13 )d 1 {<j) s d(j) s ) 

+(63 + ^6 + 2&H + 613 + 2617)9 2 (<f> s d 2 (f > s ) + (64 + 67 + 612 + 2615)9 2 {d(j ) s ) 2 
+' 5 < /'i '|(^2 + 2&4)/> s 9 1 4 > s + (65 + 267 + 615)9 1 ((f> 2 ) + bgd 1 d(j) s d 1 <f> s 

+(69 + 26 i 2 ) 9 - 2 (<^ 9 <^) + 6 i 3 9 - 2 ( 9 2 (/» s 9 -V s )| 

+d 2 cj)i{b 3 (l) s d~ 2 (t) s + 6 6 9 _ 1 (^ s 9 _ Vs) + 6 n 9 _2 (^ 2 ) + 6 i 3 9 ~ 3 (<)> s 9</> s )1 - [E.V.]. ( 47 ) 


Here the subscript ‘SCV’ indicates that this is a super cosmic variance contribution to the quadratic term and we have 
used the conditions from Eq.(41) to remove some of the original 18 6^. By design, the largest cosmic variance terms 
(those terms proportional (f>i ) regenerate the four quadratic terms which span the local, equilateral, and orthogonal 
bispectra types. The contributions proportional to d(j>i or 9 2 />; are subleading, and likely to be unobservably small 
(although we will comment further on those terms in the next section). Focussing on the leading contributions, we 
can more clearly make contact with our discussion of bispectral shapes (Eq. ([9]) ) by writing 


$ 


SCV 


gNhfii 

/nl 


{[^i^s+^29 l ((f> s d(f> s ) + A 3 d 2 {(j) s d 2 (j) s ) + A^d 2 (9(/> s ) 2 ] + • • ■ — [E.V.]| 


(48) 


(s) 

There is also a term in Eq.(44), 4>3 , with a single <j> a : the cubic terms also contribute to the linear term in biased 
sub-volumes. Including the shift to the linear term from <f> 2 (see Eq. (32)), the induced linear term is 


$? CV = /nl 


+ffNL 


</>s{(2ai + a 2 + a 3 )(fn + (a 2 - 2 a 3 )d(f>i - + a 3 d 2 <f)i—^ 

0s'|[36i + 6 2 + 63 + 65 + &6 + 26 8 + 6n + 2 bn ]( j ) 2 + [6 2 + 69 + 613)9 1 ( cf > id ( t > i ) + 

b 3 d~ 2 {(j)id 2 (l)i) + 6 4 9' 2 (90i) 2 


+9 4 > s < [6 2 + 264 + 265 + 267 + 69 + 26 i 2 + 2615] ( rf ’ idrfu ) + 6 6 9 ( 0 z 9 <fii) + b^d (d(f>i) + [6g + 6i5]9(/> ; ) 


+ 9 2 <M [63 + 6 6 + 26 n + b l3 \{(j)id 2 (j>i) + [69 + b 12 ](d4n) 2 + b 17 d 2 ((/)if 


^ [E.V.] 


9 b 13 (d-(j)id(j)i) 


(49) 


All the terms inside the first set of curly brackets in the <?nl term generate SCV of the same order (that is, 
((^?) 2 ) ~ ((9 _1 (</>z9</>z)) 2 ), etc), while the terms on the fourth line and below give contributions that are significantly 
smaller. Furthermore, since the individual terms inside each set of (?nl curly brackets have the same dimension, there 
is no way to distinguish them. The only relevant feature is the degree of scale-dependence of the prefactor. 

Finally, notice that the short-wavelength field contains terms where, for example, p±, p 2 , p 3 > k* but |pi +P 2 I < &*• 
This momentum range allows us to consider the limit of correlation functions where momentum modes accessible in 
a sub-volume sum up to a long wavelength mode. The information in that limit is about the variation of lower order 
correlations over spatial distances the size of the long wavelength mode mm 


IV. FAMILIES OF CUBIC TERMS FROM THE SUPER COSMIC VARIANCE POINT OF VIEW 


In this section, we point out various interesting special cases of the general formulae from the previous section. 
In each of the subsections below we impose additional constraints (beyond the loop correction considerations) that 
restrict the super cosmic variance contributions to the bispectrum and power spectrum. That is, we impose constraints 
on the various combinations of the 6^ parameters in Equations (47) and (|49|) . 
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A. General results 


Before considering the results of the previous section organized by the squeezed limit behavior of the bispectrum, 
we first note some general features of the possible super cosmic variance consequences from our cubic terms. 


• We find cubic terms that shift the bispectrum at order (j>i but that give no shift to the power spectrum at order 
(j> a (the terms in Eq.(W9j) proportional to </>/, <9 _ 1 (<(>z<9/>z), d~ 2 (c/)id 2 4>i) and d~ 2 {d(j)i ) 2 individually vanish). This 
set of cubic terms only"generates $ 2 CV where the {Ai} are linear combinations of { 1 , 0 , — 2 , 2 } and { 1 , — 2 , 0 , 0 }. 


• We find cubic terms that shift the power spectrum at order <p s but that generate no significant shift to the 
bispectrum (all terms proportional to (f>e in the induced quadratic term individually vanish). The $f cv generated 
by any terms in this set are indistinguishable at order cj) s . Some can, in principle, be distinguished by the relative 
amplitude of the sub-leading terms (proportional to d~ 1 cj) s and d~ 2 (j) s ). 


We find cubic terms that generate neither a shift proportional to cj)(, in the induced quadratic term, Eq.(47), nor 
a shift proportional to (f> s in the induced linear term, Eq.(49). From an observational point of view, these cubic 
terms generate no super cosmic variance of an interesting size. One of these induces no sub-leading shift to the 
power spectrum and at most a d 2 <fii shift to the bispectrum and so is a candidate for a trispectrum template 
consistent with single clock inflation. The kernel and trispectra for this case can be found in Appendix |D| 


To summarize, our set of cubic terms demonstrates that it is possible to find examples with significant cosmic 
variance in the bispectrum but not the power spectrum, in the power spectrum but not the bispectrum. As may be 
expected, our set also contains terms that generate no interesting cosmic variance. Furthermore, there is, in general, 
a degeneracy in the induced lower order terms: Multiple cubic terms (eg, terms that give different trispectra) can 
generate indistinguishable 4>® cv , $ 2 CV . 


B. Cubic terms that induce equilateral type bispectra in biased sub-volumes 


If we impose the conditions 4Ai + 2A-± + 2A 3 = 0, and 2A 2 — 4 A 3 = 0, the SCV induced quadratic term will generate 
a bispectrum of the equilateral shape (the condition A 3 + A 4 = 0 is already enforced by the loop constraints on $ 3 ). 
Most cubic terms satisfying those constraints shift the power spectrum at leading order. However, there are cases 
that generate an equilateral bispectrum but no leading order shift to the power spectrum (no term proportional to 
4> s in Eq.([49|). This case may, in principle, be distinguished from a purely equilateral bispectrum by a sub-leading 
term (proportional to <9</>j), although in practice that contribution is quite small. There is also a cubic term whose 
SCV contribution gives an equilateral bispectrum at leading order and no contributions to the bispectrum of size dfa 
(although this solution does affect the power spectrum at order (f > s ). 

We can define an equilateral family of statistics by those whose SCV contributions induce only an equilateral 
shape bispectrum at leading order and shift the power spectrum by at most terms proportional to d~ 2 (f> s (the same 
effect that an equilateral bispectrum gives). Note that although this family is distinct from the set of correlations 
with no SCV at all, it is not observationally distinct if only the power spectrum and the shape of the bispectrum have 
been measured. A set of non-zero { 6 j} for the cubic term that induces the equilateral bispectrum but no significant 
shift to the power spectrum is 


Equilateral family cubic term : b\ = —5/3 , b 2 = 2 , 65 = 3 , 69 = —4 , = —2 , &13 = 2. (50) 


The field observed in biased sub-volumes in this scenario (assuming only an equilateral bispectrum in the mean 
statistics) takes the form 

4> obs (k) = 0(k) (51) 

/ /»equil | equil / \ /. 

+ /nL 2 !( 9 ^3 L —— J d 3 p 1 d 3 p 2 [<Kpi)<Mp2 ) - (<HPi)<HP2)>] N 2 eqUll (Pi; p 2 , k)£ 3 (k - Pi - p 2 ) 


equil 3 » 

/a, 


0(pi)0(p 2 )</(p 3 ) - ^2 <HPiK0(Pj)<MPfc)) 


i= 1 


-^ 3 quil (Pij P 2 , P3; k)<5 3 (k — ^2 Pt) 


where N ^ qml is Eq.(|lo| with { 0 ^ = {—3,4,2,—2} (see Eq.(20)) and N^ qml is Eq.(A5) with the set of {&;} from 
Eq.(|50|). Figure [ 2 ] shows an example quantification of the qualitatively important result for this scenario: there is a 
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range of values of i n the mean or parent statistics that are consistent with a particular value of J^ ull ' obs j n 

sub-volumes. Note that the parent volume can have an input < j! ml of either sign. Appendix D contains more details 
of some of the interesting cubic terms that induce an equilateral bispectrum in biased sub-volumes. 


pi/nT') 



FIG. 2: Assuming that a single sub-volume measurement of f ^ 11 = 10, this plot shows the probability distribution for the 
value of that could be present in the (mean statistics of the) parent volume. Both curves assume the parent volume has 
a power spectrum consistent with Planck satellite observations and g)^ 1 ' = 5 x 10 3 (where g^ 1 ’ is defined as the amplitude 
of the non-local cubic term specified in Eq.(50l). If there are no further SCV effects from higher order terms, any sub-volume 
will have the same power spectrum and g ^) 1 as the parent volume. The solid (blue) curve considers a parent volume larger 
than the sub-volume by a factor of 50 extra e-folds. The dotted (black) curve show a parent volume larger by a factor of 100 
extra e-folds. For the best fit Planck power spectrum and N e = 50, 100, the RMS variance of the long wavelength modes is 
y/{(jyf) ss 0.0008 and 0.0022 respectively. 


C. Cubic terms that induce local type bispectra in biased sub-volumes 

If we impose the conditions A 2 = A 3 = A 4 = 0, the SCV induced quadratic term will generate a local shape 
bispectrum. All such cubic terms also shift the power spectrum at order <p s and most also give sub-leading shifts. 
If we define the local family as the set <h n whose SCV contributions induce only a local bispectrum and only shift 
the power spectrum as cf> s (no sub-leading terms, to match the behavior of the local bispectrum), our ansatz contains 
several cubic terms in the family. If we are more restrictive and insist that the induced bispectrum has no piece 
proportional to either d(f>i or d 2 (f>i in the bispectrum, the only solution is bi = 0 except for i = 1 (the local g^L term). 
Super cosmic variance from a local cubic term and beyond was previously considered in [2J I2'6] . 


D. Cubic terms that induce orthogonal type bispectra in biased sub-volumes 

If we impose the conditions 2Ai + A 2 + A 3 = 0, the SCV induced quadratic term will generate a bispectrum with 
a squeezed limit that diverges as 1/kf (the condition A 3 + A 4 = 0 is already enforced by the loop constraints on 
<I> 3 ). All such cubic terms in our set generate no shift to the power spectrum at leading order. Notice that this does 
not insist on the orthogonal template, only the squeezed limit behavior. However, it is possible to have the exact 
orthogonal template bispectrum induced with or without also inducing a shift to the power spectrum. As in the 
local and equilateral case, we can define the orthogonal family as the set of terms in [gf>(x)] = </>(x) + /nl^M^x)] + 
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5 nl < £ > 3 [ ( /K x )] + • • • that has no stronger than <p s /k SCV in the power spectrum and maintains a bispectrum with only 
1 /kf divergence. 


V. DISCUSSION 

In this work we have considered non-Gaussian fields built from non-local real space expressions with up to cubic 
dependence on a Gaussian field. This ansatz, Eq. 0 . can generate a variety of tree-level bispectra and trispectra, 
including shapes consistent with either single-clock or multi-field inflationary models. By examining the ways in 
which statistics in spatial sub-volumes can differ from the mean statistics of the parent volume we draw the following 
important conclusions: 

• Terms at order n in a non-local expression for a non-Gaussian field $(x) can generically shift the amplitudes 
and momentum dependence of the power spectrum and correlations up to order ($(fci)...$(fc n )) in biased sub¬ 
volumes, but need not shift all of them. 

• More specifically, super cosmic variance effects can induce bispectra in biased subvolumes to match any of the 
equilateral, orthogonal, or local templates. The amplitude and sign of the induced bispectrum depends in a 
degenerate way on the background over- or under-density of the sub-volume as well as the amplitude and sign 
of the mean trispectrum. Depending on the specific form of the trispectrum, the power spectrum may or may 
not display a dependence on background density within the sub-volume. 

Mathematically, our results are very reasonable: non-Gaussian statistics at each order are independent, so of course 
measuring the 3-point function alone constrains neither the shape nor the amplitude of higher order correlations. 
Furthermore, different n-point correlations can induce indistinguishable contributions to lower order correlations in 
biased sub-volumes. When neither the bias nor the higher order correlations can be measured, there is a large 
degeneracy of correlation functions (and so inflation models) that can be consistent with just a few measurements. 
Specifically, we have demonstrated that a limited set of measurements of the primordial density correlations (e.g., 
a detection of /^L al = 0 and /n't )" 1 7 ^ 0 ) could be consistent with single-clock inflation or a multi-field scenario. 
Measurement of a purely equilateral bispectrum does not in and of itself imply that there can be no super cosmic 
variance at work. To rule out super cosmic variance (at least, to rule it out of observational relevance) we must 
also constrain any correlation between statistics measured within smaller regions of our own Hubble volume and the 
background density of those regions. 

We have hardly covered the space of non-Gaussian statistics: our set of examples was chosen to be a minimal 
set that allows us to explore the possible implications of cosmic variance from super horizon modes for non-local, 
scale-independent bispectra 3 . We see no obstruction to extending our results to find quartic terms in the non-local 
expansion that could bias, for example, the trispectrum but not the bispectrum. One might complain that super 
cosmic variance that evades detection in the simplest measurements (e.g., in the power spectrum) but biases higher 
order terms (e.g., the trispectrum) would require too much of a conspiracy to be realistic. Addressing this discomfort 
systematically requires a measure on the space of super cosmic variance effects from multi-field inflation models. It 
would be interesting to find an inflation model that contains trispectra of the sort we have found or to prove that no 
model can generate them, nor their higher order generalizations. 
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3 There are many other terms one might wish to include to study further aspects of this problem. At the level of the quadratic terms, for 
example, we have not been sufficiently general to cover the other standard template for single-clock inflation (the ‘orthogonal’ template 
in the appendix of reference [23]). To expand the model to include this template would require allowing a.i < — 2. At cubic order, 
we eliminated some terms since they did not induce one of our desired quadratic terms at leading order. Keeping those terms would 
add more parameters that would relax some of the restrictions on the bi from removing sub-leading terms (and badly behaved loop 
contributions). 
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Appendix A: Cubic terms and the cubic kernel 


The functional $ 3 [</>(x)] built from the 18 terms derived in Section III A 


is 


$ 3 (x) = bicj) 3 + b 2 <j)d~ 1 (<j)d(f>) + b 3 4>d~ 2 ((/)d 2 cj)) + 6 4 ^<9 -2 ((<9<^>) 2 ) + 6 5 9 _1 (^ 2 9^) + b^ 1 (c/)d 2 cf>)) 

+b7d~ 1 (<^<9 _1 ((dc/)) 2 )) + &s<9 -1 {(j)d{(j) 2 )) + b 3 d~ 3 (dcjid -1 {(j>dcj))) + 6i 0 <9 -1 (dcj)d~ 2 (cj)d 2 (/>)) 

+bnd~ 2 (cf 2 d 2 (j)) + 6i2<9~ 2 (</>(<9(/>) 2 ) + bi 3 d~ 2 (d 2 (jid- 1 {(jidcj))) + bud~ 2 [d 2 4>d^ 2 (4>d 2 cj))) 

+bi 5 d~ 2 (d(/)d{(l) 2 )) + b w d^ 2 ((j)d{(j)d(j))) + &i 7 5 _2 (05 2 (^ 2 )) + b 13 d~ 2 (dc^d -1 {cj>d 2 $)) . (Al) 

To compute the kernel associated with 4> 3 (x), we need to take the Fourier transform of each term. As an example, 
consider the Fourier transform of the 13th term: 




= d 


l 


k 2 

d 3 pi 

(2tt) 3 


%p50(pi)e‘"' 


J (271-) 

d 3 p 2 

(2tt) 3 


d 3 q 1 

(27t) 3 q . 


d 3 p 2 

(2tt) 3 


d 3 p 3 S 3 (q - P2 - P3)^(P2)P3<A(P3) 


-ik : 


d 3 p 3 J d 3 q <5 3 (k - pi - q)<5 3 (q - p 2 - P 3 )^|^</>(Pi)</ > (P 2 )</>(P 3 ) 


(A2) 


We can use a Dirac delta function to integrate over any of the four momenta (pi,P 2 ,P 3 or q) but to make the 
calculations more transparent we will treat these possibilities symmetrically and write 


PiPl 


+ 


-11 

' d 3 p 1 
(2tt)3 . 

f d 3 p 2 

1 (27T)3„ 

j d 3 p 3 S 3 {k 

P2P 2 3 

-I- 

PiPl 

PlP3 


P 2 lP2 


p\P3 


fc 2 |Pl+P 2 | k 2 IP 1 +P 2 I & 2 |pi+p 3 | fc 2 |pi+p 3 | fc 2 |p 2 +p 3 | k 2 |p 2 T P31 


(A3) 


Now we can read off the contribution to the cubic kernel from the 13th term: 


^ r 3 13th (Pl,P2,P3,fc) = ' 


P2P 2 3 


PiPl 


p\p3 


+ 


P 2 lP2 


p\P3 


k 2 IP 1 +P 2 I k 2 |pi T P 2 1 A 2 |pi+p 3 | fc 2 |pi+p 3 | k 2 1P 2 T P 3 1 k 2 |p 2 +p 3 


(A4) 


4 We note that three terms - d~ 1 (d(f)d~ 2 (dcj)) 2 ), d~ 2 (d 2 4>d~ 2 (<9</>) 2 )), d~ 2 (d(/)d~ 1 (d(j)) 2 ) - have been neglected because they do not 
have nonzero 0±, 02, or 04 and thus do not regenerate the original four quadratic terms, but only terms suppressed by d(f>i. Allowing 
these terms would introduce new parametric freedom among the 18 terms above, allowing some to be nonzero after the imposition of 
convergence on loop corrections. 
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Applying this procedure to each term gives the full cubic kernel iV^pi, p 2 , P 3 , k): 


,|Pi +P 2 I 
P 1 P 2 


^3(Pl,P2,P3 ,k) 

ru L h I Pi +P2 

= 661 + 02 
+ 25 4 
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iPl +P2| 


k 2 


IPl +P3| 
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+ 


k 2 


IP 2 + Psl 
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k 2 


(pj+pj)P3 (Pj +pj)P2 (Pl+PDPI 
k 2 1P 1 + P 2 1 k 2 1 p 1 + P 3 1 k 2 1P 2 + P 3 1 


(A5) 


Appendix B: Example of finding the quadratic terms induced by a cubic term in biased sub-volumes 


We again use the 13th term, d~ 2 (d 2 (jid -1 (4>d(j))), to illustrate how to read off the induced quadratic terms from 
the limit of a cubic term when one momentum is much smaller than the other two. Choosing any one of the three 
momenta (e.g., p%) as the long wavelength mode ki and the other two (pi,P 2 ) as the short wavelength modes k s , and 
using the fact that the kernel is symmetric in the pi: 


d 2 (d 2 4>d (4>d4>)) (k) 

d 3 P: 


3 x ^ 

3! J a 


(27r) 


3 </>(P3) 


f 

Ik, 




Pi 


i d 3 p 1 

(2tt) 3 

P 2 


d 3 p 2 d 3 (k - pi - p 2 )<?i(pi)0(p2) 


<k. 


Pl 


pi 

k 3 


El 

k 3 


k 2 p 2 k 2 pi 

\ / 

= 1 4nd~ 2 ((/) s d 2 (l) s ) + d(jnd~ 2 (d 2 <j) s d- 1 (t> s ) + <9 2 <(>z<9 _3 (</i s <9</> s ) 


(Bl) 


where A is the largest scale in the problem (eg, corresponding to the mode that first exited during inflation) and k * 
is the scale that defines the size of the sub-volume. The corresponding procedure in real space is 


d “( d 2 cj)d (<l>d(j))) 


d~ 2 (d 2 (l>id~ 1 ((j) s d(j) a )) + d~ 2 [d 2 (j) s d~ 1 ((j)id(j) s )) + < 9 ~ 2 (d 2 0 s < 9 _1 (<(> s < 9 (/>;)) 
<9 2 <^<9 _2 (d _1 (<(> s <9^ s )) + (t)id^ 2 (di 2 (j) s d~ 1 (d(j) s )) + <9<^<9~ 2 (<9 2 ^ s <9 _1 (</> s )) 
cj)id~ 2 ((j) s d 2 (l) s ) + d4>id~ 2 (d 2 (j) s d~ 1 (l) s ) + d 2 ct>id~ 3 ((j) s d(l) s ) 


(B2) 


Note that to go from the first line to the second, all the long modes and d 2 (\n are treated as constants. 
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Appendix C: The contribution to the power spectrum from the cubic held <f >3 


In this Appendix we derive the constraints on the terms in $3 from their contribution to the power spectrum: 
J$ 3 (k')$ 3 (k)} = SP^(k) 5 3 (k' + k) 


/ gNL 


d 3 pg j-j f d 3 qr, 


n 


V 3! ) l\J (27T)3 11 J (2ny 


;^ r 3(Pi,P2,P3,k)A 3 (qi,q 2 ,q 3 ,k , )(5 3 (k' - ^ q m )^ 3 (k - ^ p*) 


f=i 


</>(pi)</>(p2)^>(p 3) - ^(p*)(</ > (Pj)</ > (Pfe)) </>(qi)</ , (q2)^(q3) - 


i =1 


i '=1 
k'+j'+i 


~ ffNL 


r d 3 pi d 3 p 2 
(27t) 3 (27t) 3 


rf 3 P3^V 3 2 (Pi,P2,P3,k)P(pi)P(p 2 )-P(p3)^ 3 (k-pi p 2 -P 3 ) ■ 


(Cl) 


In order to obtain the constraints on the terms in N 3 we need to find those that give divergent integrals in the 
expression above. The divergences are in the UV, where some of the internal momenta are large, and there are several 
configurations to consider: 

For k < |pi + p 2 |,Pi;fc >p 2 : 


614 = 0 
bio + bis = 0 

f> 3 + &4 + bo + 67 + 2fon + 6 12 + bi 3 + foi 4 + 2615 + bie + 2&17 + 6ig = 0 
For k < jpi + p 2 |,Pi,P 2 : 


bi 2 + bi3 + bie = 0 

bn + bi 4 + 645 + bn + bn = 0 


For k < pi,P 2 ] k > |pi + p 2 |: 


bie = 0 

(bn + &15 + bn) + (&12 + ^13 + bie) + {bn + 644 + 645 + 647 + bis) = 0 

be + b’j + 6 18 = 0 

(265 + be + 2bs + bg + bio + & 12 ) + {bn + ^43 + bie) = 0 
We have verified that other loop corrections (e.g., to the bispectrum) do not give additional constraints. 


(C2) 


(C3) 


(C4) 


Appendix D: Examples of special cubic kernels and the corresponding trispectra 


A cubic kernel that gives no cosmic variance contributions stronger than d 2 4>i to the power spectrum or bispectrum 
has {bi = 1, bn = 3, 647 = —3}. This is a candidate for a single-clock inflation trispectrum: 


(Dl) 


k 2 + kg + k 2 

— ft 

" |k 3 + k 2 | 2 + |k : + k 3 | 2 + |k 2 + k 3 | 2 " 

k 4 

U 

k 2 


single clock candidate : A 3 (k l7 k 2 , k 3 , k^) =6 + 6 

Recall that the trispectrum is related to the cubic kernel by 

Pt(ki,k 2 ,k 3 ,k 4 ) = 5NL P 0 (fci)P^(fc 2 )P 0 (/c 3 )A 3 (k 1 ,k 2 ,k 3 ,k 4 ) + eye. . 


(D2) 


For the kernel above, the trispectrum has no factors of k, + k ; in the denominator and lim/ Ci ^o Pt> = 0. 

A cubic kernel that induces air equilateral bispectrum and a leading order (oc (f> s ) SCV shift to the power spectrum 
is {64 = —3, & 2 = 4, 6 3 = 2 ,& 4 = —2}. We start writing the non-Gaussian field up to cubic order with implementing 
our parameter set for the cubic terms 


Induce equilateral bispectrum and shift PS: 
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ki + k 2 k\ + k 3 


|ki +k 2 


|ki +k 3 


k 2 + k 3 
|k 2 + k 3 | 


A^ 3 (ki, k 2 , k 3 , k 4 ) = —18 
+2 
-4 

One can see directly from the trispectrum that this term induces the equilateral bispectrum: 

[.2 


r ki + q 

k\ + kl 

k 2 + fc| 

Ukr+k 2 p 

ki + k 3 2 

|k 2 + k 3 1 2 _ 

k ± k 2 

kik 3 

k 2 k 3 

[|k!+k 2 |2 

|kr+k 3 p 

|k 2 + k 3 | 2 _ 


ki lim T$(ki,k 2 ,k 3 ,k 4 ) = P{k 1 )P(k 2 ) 

& 4—>-0 


-6 


ki+ k 2 k\ + k'2 k\k 2 
- 4-:- i ^-—- — 4 




k 3 J 


perm. 


(D3) 


(D4) 


The cubic term with {bi = —5/3, b 2 =2, b 5 = 3, bg = —4, b\ 2 = —2, b 43 = 2} also induces the equilateral bispectrum 
but does not shift the power spectrum at leading order (oc <j> s ). This kernel is 


jV 3 (k-|, k 2 , k 3 , k 4 ) — —10 + 2 
-4 

-4 


k\ + k 2 k\ + k 3 k 2 + k 3 
ki+k 2 | |ki+k 3 | |k 2 +k 3 

k 3 (ki+k 2 ) k 2 (ki + k 3 ) ki(k 2 + k 3 ) 


+ 6 


k\ + k 2 + k 3 


ki |ki 


fc 4 |ki + k 3 | & 4 jk 2 + k 3 | 


kik 2 + k\k 3 + k 2 k 3 
k\ 


+ 2 


fc|(fcl + fc 2 ) fc|(fci+fc 3 ) fci(/C2+fc 3 ) 


_fc||ki+k 2 | fc||k 4 


+ 


fc 4 |k 2 


In this case the trispectrum again satisfies Eq.(D4) but in addition 

No shift to power spectrum : kfk 2 lim Tj>(ki, k 2 , k 3 , k 4 ) = 0 . 


k 1 ,k2~>-0 


(D5) 


(D6) 


In the above examples, we have not normalized the coefficients bi. However, for some template definitions it may 
be useful to require JT bi = 1, which can be easily done. 
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